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I. INTRODUCTION 

o : 

In the investigation of nonlinear integrable field and lattice systems in (l+l)-dimensions an important role is played 
by the algebra of symmetries. For these systems for which additionally exists a so called recursion operator with 
vanishing Nijenhuis torsion, the algebra of symmetries takes the form of the hereditary algebra (centerless Virasoro 
algebra). Nevertheless, the construction of the recursion operator of a given integrable system is not an easy task. 
Several methods were proposed of greater or lower generality, see for example refs. jjj, H, ||, jl. Perhaps the best 
known is the one based on the implectic-symplectic factorization of the recursion operator, i.e. on a bi-Hamiltonian 
property of the system cosidered. But then the problem is shifted to the construction of two Poisson structures 
which is a complicated problem and on a general level requires advances tools ||. Moreover, it fails in the case of 
HH ' non-Hamiltonian reductions. 

C/3 \ Recently, amasingly simple and general approach to the construction of recursion operators was proposed by Sokolov 
et al. The method is based on Lax representation of a given hierarchy of integrable systems and allows a 

* ^ | construction of a recursion operators through really elementary calculation. 

In the paper presented we develop the idea of Sokolov, applying it to different classes of Lax chains for field and 
l— "™ lattice systems. Moreover, we construct from Lax representation in a systematic way conformal symmetries, beeing 
I . complementary ingredients in the construction of hereditary algebras. 

> : 

<N ■ 

Cn ' II- PRELIMINARIES 

o : 

Let us consider the following scalar Lax operators (together with their admissible reductions): 
(i) L = d» + u N ^- 2 + u N - 3 d?- 3 + ... + uo, 
"5 ■ (ii) L = flf + UN^d?- 1 + u N ^- 2 + ... + u + d- l u- x , 
"3 ("i) L = u%d% + un^- 1 + u N - 2 d?- 2 + ... +u Q + d~ l U-i + d x 2 u^ 2 , 
for field systems and 
(iv) L = £ N+a + u N+a ^ 1 £ N+a - 1 + ... + u Q + u^£- 1 + ... + u a £ a , 



> 

X 



for lattice systems, where — N < a < — 1, N + a > 1. 
Here and further on we use the following notation for differential and shift operators 

d x a(x) = a x + ad x , D x a(x) = a x , 
£a(n) = a(n + 1)£, Ea(n) = a(n + 1). 

The related Lax equations are 
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(<) L tq 
(m) it, 



L« > ,L 
>i,L 
L% > 2 ,L 
L»fc > , L 



9=1:2, 



S-oo an d i ™ >r = X)r u »^:c- Here «i are some concrete 



where L « is a pseudodifferential series of the form L~ 
functions depending on the coefficients of L. 

Fixing N (or N and a), each case gives an infinite hierarchy of commuting flows. The case (i) was first considered by 
Gelfand and Dikii [Q, cases (ii), (Hi) by Kupershmidt ||] and the case (iv) by Blaszak and Marciniak respectively. 
For arbitrary flow from (i) — (iv) we are going to construct a respective algebra of symmetries. What we really need 
for this construction are two invariant objects, i.e. the related conformal (scaling) vector field a and a recursion 
operator <f> 



Ck<J 



da 
Of 



= 0, C K <j> = Q, 



(1) 



where if is a vector field of a given flow and C means a Lie derivative. 
Lemma 1 
If 

C a K ~ pK, C a 4> — a 4>, P,ct = const, 



and 



C^rCf) = 4>c T 4>, 

for an arbitrary vector field r, then K n and a n form a hereditary algebra 

[K n , K m ] = 0, [cr„, K m ] = (p + am)K n+m , [<r n , a rn ] = a(m - n)a n+r 



(2) 



(3) 



(4) 



The proof the reader can find for example in ref. |Q . Here we show how to extract scaling symmetry and a recursion 
operator directly from Lax equations. 



III. SCALING PROPERTIES OF LAX EQUATIONS 

Let us start from scaling symmetries for field systems. Lax equations (i) — (Hi) are homogenous with respect to 
the following scaling 



d x -> e e d x , u N -i -> e^UN-i, d tq -> e qs d tq , L -> e Ne L. 
The conformal point transformation takes the form 



(5) 



ttJV-1 



Uo 
U-i 

v u_ 2 y 



/ e 0e Mw (e £ z,e 9 %) \ 
e le WAr_i(e e x, e qe t q ) 

e Ne u (e e x,e qE t q ) 



(6) 



so the related infinitesimal generator, i.e. conformal symmetry is 
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I \ 

V (7V + 2)u_2 y 



Un-X 
Un-2 

v ui 2 y 



9*9 



/ U N \ 

V «i a y 



All admissible reductions preserve this scaling property. 

Lax representation (iv) for lattice systems is homogenous with respect to the following scaling 

S^e £ S, u N+a ^ ^ e l£ u N+a ^, d tq ^ e qe d tq , L^e {N+a) L. 
The conformal point transformation and conformal symmetry are 

/ u N+a -i \ ( e l£ UAr +Q _i(n, e q H q ) \ 



U 



(7) 



(8) 



I UN+a-l \ 

{N + a)u 



UN+a-2 



\ Nu a J \ 



(9) 



IV. ALGORYTHMIC CONSTRUCTION OF RECURSION OPERATORS 

Developing the idea of Sokolov et al. || we construct recursion operators directly from Lax hierarchy 

L tq = [A q ,L}. (10) 

The only necessary information is the explicit form of L operator and the ansatz A q = A q P + R relating A operators 
with different q, without any specific information of the coefficients of A. Here, P is some operator that commutes 
with L of the form P = P(L) and R is the remainder, hence 

L tj = L ti P(L) + [R,L]. (11) 



Remark 

In fact, formulas similar to ( pi] ) appeared for the first time in 80th, derived on more general level in the frame of 
r-matrix formalism (see for example refs. poj , jnj, O], [^3|). Nevertheless, they contain r-matrices and involve more 
calculations. The advantage of presented method is its simplicity, although the roots in r-matrix theory are evident. 

The case (i), i.e. the Gelfand-Dikii case, was considered in details in ref. |J, here we concentrate on the cases (ii) 
and (Hi) for field systems and on the lattice scalar Lax equations (iv). 

Of course, in order to apply recursion operators, found by the method presented, to the construction of hereditary 
algebra of symmetries, it is necessary to verify their hereditary property (|J). 
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A. Field systems (ii) 



(N + l)-field Lax operator and the Lax hierarchy are 



u 



d x 1 U-%, 



L tm = [A m ,L], A m = (L$)> 1 , m=l,2,... . 

Let us express the A m+ n operator through A m , L and some remainder R m 

A m+N = {L^L)>i = (L® > X L + Lt? <x L)>i 
= > X L - (L® > X L) Q + (L% <ii)>i 
= A m L + R m . 



(12) 



(13) 



Analysing the highest and lowest order terms of R m we conclude that the remainder is a purely differential operator 

R m = a m + b m d x + ... + 7 m 3f, (14) 

hence 

= L tm L + [R m ,L]. (15) 



The right hand side of eq. ( |l5| ) is the pseudodifferential operator: RHS — v 2 N-id x N ^ 1 + ...+v +d~ 1 v-i +w_ 2 5~ 1 it-i. 
From Vk = for k = 2N — 1, N, — 2, we determine R m coefficients in terms of the coefficients of operators L and 
Lt m . Comparing the remaining coefficients of both sides of operator equation (fiq) we get the recurrence formula 



(16) 



/ UiV-1 ^ 




( u N - 1 \ 


V «-! J 


-* 

tm + N 





There are iV commuting chains of vector fields generated by the recursion operator 

K r>n = 4> n K r , r = l,...,N, n = 0,l,2,... . 



(17) 



Let us pass to the admissible constraints. There are few of them 14 1, ||. The first constraint is of the form u_i = 
and is non-Hamiltonian. It means that we cannot apply a bi-Hamiltonian method for construction a recursion operator 
and so the method presented becomes even more important. The second constraint takes the form u_i = ua = and 
preserves the Hamiltonian structure. The last constraint is the so called Kupershmidt reduction M 



L = (-l^d^tfda;, q-odd => u N -i = 0, 



(18) 



where (ad£)^ = (— l)"9"a. Notice that under this constraint half of equations from the hierarchy ( [To|) disappear and 
only these with odd q remain. Two cases have to be considered. The first one with an even N preserves the recurrence 
formula (|l5|) as odd m and even N give odd m + N. The second case of odd iV is more complex as m + N is even 
and so it is excluded from the hierarchy. In this case we must take 



A m+ 2N = (LL^L)>! = {LL« >ii)>i + {LL« <x L)>i 

= LL^ > X L - (LL^ >iL) - (LLt* >iL)_i + (LL% <i£)>i 
= LA m L + R m , 



(19) 



where R m = a m d^ N 1 + ...+d x lr ) m . In derivation of the formula ( |l9| ) we applied the identity d x 1 a x d x 1 = ad x 1 — d x 1 a 
and the fact that Ujv-i — 0- Hence, the new recurrence formula is 



L t m+2N = LIj tm L + [R m , L], to - odd. 
Example 1. Consider Lax operator of the Kaup-Broer system 

L = d x + u + d~ 1 v, 



(20) 
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with the remainder R m = a m d x + b m , hence 

L>t m+1 = L tm L + [R m ,L] 
t 

u tm+1 + d x x v tm+1 = [u tm - (a m ) x ]d x + [uu tm + v tm + a m u x - [b m ) x ] 
+d x 1 [uv tm - (v tm ) x - v(D~ 1 v tm ) + (a m v) x ] 
+ [u tm + {D- X v tm ) + b m ]d x 1 v. 

Finally 

[ut m ~ (a m )x] = => a m = D^ 1 u tm , 

[ u t m + {D x X v tm ) + b m ] = => b m = -u tm - D~ 1 v tm 
and the recursion formula ([l6]) takes the form 

u\ ( D x -\- D x uD~ l 2 \ ( u 



v I \ v + D x vD T 1 -D, x + u ) \ v , 

Imposing the constraint w = 0we get the Burgers hierarchy with 

L = d x + u, a m — > D^u^, b m — > —u tm 
and the recursion operator of the form 

c/) = D x +u + u x D- x = D x + D x uD~ x . 
The Kupershmidt reduction with odd N = 1 gives the nonstandard KdV Lax representation 

L = d x + d x 1 v 

with Lax recursion 

Lt m+2 = LL tm L + [Rm, L], rn — odd, 

where 

Rra dm&x ~t~ b rn ~t- d x C m , 

a m = D~ 1 v tm , b m = 0, Cm = -D x v tm - vD~ 1 v tm 
and related vector field recursion 

«W, = (D 2 x +4v + 2v x D- 1 )v tm . 
Example 2. Consider a three field Lax operator 

L = d x 4- ud x + v + d x 1 w. 

The appropriate remainder takes the form 

Rm — @"m& x ~t~ b m d x -\- C m , 

where 

a m = 7^D~ 1 u tm , b m = ^uD~ 1 u tm - jD~ 1 uu tm + ^D~ l v tm , 
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Cm = ~V tm - D x 1 W t „ 

So the recurrence formula (n6h is 



u 
v 
w 



A \D x + \D x uD x x 3 
B C 2u 

D jw + ^WxD^ 1 Dl — D x u + vJ \w 



where 



A = \d\ + D x vD x x - -DxuD^u, 

3 2 1 1^,1, 

B = ~v x + -w + {w x + -v xx + -uv x )D x - -v x D x u, 

C = D x + uD x + v + -v x D x \ 

D = -jwD x + ^-uw - ^w x + ~[(uw) x - w xx \D~ x - ^w x D~ x u 



and hence two commuting hierarchies of vector fields are 

K hn = <f> n K u K x 



2v x 

K 2 , n = 4> n K 2 , K 2 =| v xx + 2w x + uv x 

-w xx + (uw) x 

The first reduction w = leads to L — d x + ud x + v, c m = —v tm (a rn , b rn are unchanged), 

\Dl + D x vD^ - \D x uDz x u \D X + \D x uD^ 

3„, i fl„. i 1„,„, \n-l 1„, n-l„, n2 , „,n i „, i 1„, n-i 



f«x + {\v xx + ±uv x )D x - jv x D x 1 u D x + uD x + v+ \v x D x 



K 1 = ( U *),K 2 = ( 2v * 

\Vx J \V XX + UV X 

The second reduction gives the MKdV hierarchy 

L = d x + ud x , 

a m = ^D x 1 u tm , b m = -uD^u^ - ^-D~ l uu tm , c m = 0, 

<\>=\d 2 x - X -D x uD x x u, K x = u x , K 2 = 0. 
Finally, the Kupershmidt reduction gives another nonstandard KdV case 

L = dl + v- id-V, 



where 



and 



Rm = a m d x + b m , a m = -D x 1 v tfn , b m = -^v tn 



= D 2 x +v + -v x D~ x 



Notice that L = <j>^ and the rescaling v — > 4i> gives the KdV case from the previous example. 



G 



B. Field systems (in) 



(N + 3)-field Lax operator and the Lax hierarchy are 

L = u N d x + v,N-idx~ + ■■■ + uo + d~ l u-i + d x 2 u- 2 , 

L tm = [A m ,L], A m = (L%)> 2 , m = l,2,.... 

Let us express the A m+ n operator through A m , L and some remainder R m 

A m +N = (L%L)> 2 = (L% > 2 L + L% <2 L)> 2 

= L% > 2 L - (L 3 ® > 2 L) - (L^ > 2 L)i + (L% <2 L)> 2 (21) 
= A m L + R m . 

Analysing the highest and lowest order terms of R m we conclude that the remainder is again a purely differential 
operator 

Rm = a m d* +1 + b m d£ + ... + 7m , (22) 

hence 

L tm+N =L tm L+[R rn ,L}. (23) 

There are several admissible constraints. The first group is the following: u- 2 = 0, u- 2 = u_i = 0, u- 2 = it_i = 
uq = 0, u_2 = it_i = uq = i*i =0. The first three are non-Hamiltonian the last one is Hamiltonian. Another 
constraint is related with the Kupershmidt reduction 

L = {-lfd- 2 tfdl (24) 



which survievs only odd terms from the hierarchy. Again an even N preserves the recurence formula (|23j) as odd m 
and even N give odd m + N, while for odd N m + N is even and so it is excluded from the hierarchy. In this case we 
must take 

A m+2 N = {LL^L)> 2 = (LL^ > 2 L)> 2 + (LL^ <2 L)> 2 

= LL^ > 2 L - (LL% > 2 L)i - ... - (LL% > 2 L)- 2 + (LL^ <2 L)> 2 

= LA rn L + R m , (25) 
where R m = a m d 2N+1 + ... + d x 2 j m . Hence, the new recurence formula is 

L t m +2N = LL t m L + [R m , L], m- odd. (26) 
Example 3. Consider 4-field Lax operator 

L = ud x + v + d~ w + d~ z. 
The remainder and the recursion operator are 

Rm Q>m & x ^w. &x ~t~ Cyn 3 

a m = u 2 D^ 1 u^ 2 u tml b m = 2D~ 3 z tm + D~ 2 w tm - u tm , 

c m = -D x 2 z tm - D x 1 w tm - v tm , 

I A u u x D~ 2 - uD- 1 2u x D x 3 - 2uD x 2 

5 v + uD x 2u + v x D x 2 2v x D x 3 +uD x 1 

(f> ~\ C w v - D x u + w x D x 2 + wD- 1 u + 2w x D x 3 + 2wD x 2 I 

\D z z x D x 2 + 2zD x 1 E 
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A = v + uD x + (u xx + 2v x )u 2 D x x u 2 , 

B = 2w + (y xx + 2w x )u 2 D~ 1 u~ 2 + 2wuu x D~ 1 u~ 2 , 

C = 'iz — 3w x — 2wD x — 2wu~ 1 u x + [2(u z) x — (u 2 w) xx }D x 1 u~ 

D = -3z x - 2zD x - 2zu~ x u x - (u 2 z) xx D~ 1 u~ 2 , 

E = v- D x u - wD x x + AzD x 2 + 2z x D x 3 . 



The hierarchy takes the form 



K n+1 = (j) n Ki, K x = _ ( " , ^, 2 



u 2 u xx + 2u 2 v x 
u 2 v xx + 2u(uw) x 
(u 2 w) xx + 2(u 2 z) x 
(u z) xx 



The first constraint z = 0, L = ud x + v + d x 1 w gives 

a m = u 2 D x x u- 2 u tm , b m = D x 2 w tm - ut m , c m = -D x l w tm - v t „ 

A u u, x ^ x 

B v + uD x 2u + v x D x 2 

-2>w x — 2wD x — 2wu~ l u x _ „_ 9 __i 

-(u 2 w) xx D^u- 2 W v-D x u + Wx D x 2 + W D x i 

(u 2 u xx + 2u 2 v x 
u 2 v xx + 2u(uw) x 
-{u 2 w) xx 

The second constraint z = w = 0, L = ud x + v gives 

a m = u 2 D x 1 u^ 2 u tm , b m = -u tm , c m = -v tn 

v + uD x + (u xx + 2v x )u 2 D~ 1 u~ 2 u 

v xx u 2 D x x u- 2 v + uD x ) ' 

u 2 u xx + 2u 2 v x 



K n+1 = <fTKi, Ki=[" " x * 2 



^ v xx 



Finally the third constraint z = w = v = 0, L = ud x gives 

a m =u 2 D~ 1 u~ 2 u tm , b m = -u tm , c m = 0, 



= uD x + u xx u 2 D x 1 u 2 , 



K n+ i = <j) n Ki, K\ = u u xx . 



C. Lattice systems (iv) 

AT-field Lax operator and the Lax hierarchy are 

L = S N+a + u N+a ^ 1 £ N+a - 1 + ... + u a £ a , a = -1, —N, N > 2, 

L tm = [A m ,L], A m = (L^)> , m=l,2,... . 
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Let us express the A m+ N+a operator through A m , L and some remainder R m 

A m+N+a = (L^L)> = {L-^h > L + <0 L)> 

= i^fe > L - (L^fe >o^)-i - ••• - {Lifts > L) a + (L^fe < i)>o 

= A m L + R m . (27) 

Analysing the highest and lowest order terms of R m we conclude that the remainder is a shift operator of the form 

R m = aj^"- 1 + ... + lm £ a , (28) 

hence 

Lt m+N+a = L tm L+[R m ,L]. (29) 
There are two possible reductions. The first one takes the form 

u a ^ 0, u a+ i = ... = u N+a -i = (30) 

L tm = [L m > ,L], m=(N + l)n, n= 0,1,2,..., 
and includes all Bogoyavlensky lattices, while the second one is 

u = it_i = ... = u a = 0, ux,...,u N+a -x ^ (31) 

L tm = L4 m ,L], A m = (i^)>i, m = 1,2, ... 

and is a discrete Gelfand-Dikii analog. 

We present the calculations on simplest well known examples of Toda and Voltera lattices. 
Example 4- Consider infinite Toda lattice N — 2, a = — 1. Hence, we have 

L = £ +p + v£~ x , R rn = a m + b m £~ 1 

and the Lax hierarchy ( p9[ ) 

Pt m+ i + vt m+l £~ X = [pt m + a m - (Ea m )}£ + [pp tm + v tm + b rn - (Eb m )] 

+ [ v Pt m + (-E^VH™ + va m - v(E^ 1 a m ) + b rn (E~ 1 p) 

-pb^ 1 + [(ET^H™ + bmiE^v) - v{E~ 1 b m )]£- 2 . 

Introducing the operator A := E — 1 and its inverse A -1 such that A _1 /(n) = X)fc=-oo /(^)? we nn d 

Pt m + a m - (Ea m ) = => a m = A~ x p tm , 

{E- l v)v tm + b m (E- x v) - v{E- l b m ) = => b m = -vA^Ev^v^ 



and then 



(E- 1 p)+v{E- 1 Ap)A- 1 Ev- 1 vil + E' 1 ) 
1 + AvA~ 1 Ev~ 1 p 

vE^ApA^Ev- 1 v(l + E- r ) 
(Ev - vE-^A^Ev- 1 p 



K - APTC K - f V ( n )\p( n ) ~P(n- 1)] 

The constraint p = leeds to the infinite Volterra system 

L = £ + v£-\ L t2m =[L 2m > ,L], to =1,2,. 







Now we find 

a (j 2m t2\ it 2m t2 , r 2m j2\ 

A-2m+2 — (-L L J> — \L > a L + L <( )L )> 

t 2m j 2 it 2m t2\ , / j 2m j2\ 
— L> >0-^ — \Li >0-k )-2 + {-Lj <0-^ J>0 

= A2, n L 2 + i?2m, 

where R^m = a m + b m £~ 2 , so the recursion chain is 

Lt 2m +2 — Lt 2m L 2 + [R2m,L] 

and hence 

a m = A _1 u t2m , b m = -v(E- 1 v)A- 1 Ev~ 1 vt 2m , 
4> = v + (E^v) + vE- 1 + Av(E- 1 v)A- 1 Ev- 1 = v(l + E~ l ){EvE - w)A" 1 w" 1 , 
K x = v(n)[v(n + 1) - v(n - 1)]. 
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